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Abstract —The response of superconducting pair-breaking de¬ 
tectors is dependent on the details of the quasiparticle distribu¬ 
tion. In Kinetic Inductance Detectors (KIDs), where both pair 
breaking and non-pair breaking photons are absorbed simul¬ 
taneously, calculating the detector response therefore requires 
knowledge of the often nonequilibrium distributions. The quasi¬ 
particle effective temperature provides a good approximation to 
these nonequilibrium distributions. We compare an analytical 
expression relating absorbed power and the quasiparticle effec¬ 
tive temperature in superconducting thin films to full solutions 
for the nonequilibrium distributions, and find good agreement 
for a range of materials, absorbed powers, photon frequencies 
and temperatures typical of KIDs. This analytical expression 
allows inclusion of nonequilibrium effects in device models 
without solving for the detailed distributions. We also show our 
calculations of the frequency dependence of the detector response 
are in agreement with recent experimental measurements of the 
response of Ta KIDs at THz frequencies. 


I. Introduction 

Kinetic Inductance Detectors (KIDs) [^]-^ rely on photons 
with energy hv > 2A, where v is the photon frequency and 
A the superconducting energy gap, to break Cooper pairs and 
thereby produce an excess population of quasiparticles. Sub¬ 
gap readout photons are also absorbed, and they too generate 
an excess population of quasiparticles, which influences the 
operating characteristics of the device i I- Understanding 
the way in which these distributions are created, and interact, 
is central to understanding the operation and performance 
limitations of KIDs. 

Once a signal photon is absorbed, the initial high-energy 
quasiparticles relax towards the superconducting energy gap 
A. For quasiparticle energies A. < E Llo vn typical 
materials, this happens primarily by emitting phonons [^. In 
a thin film, these excess phonons can either remain within 
the superconductor and break pairs themselves if they have 
sufficient phonon energy fl > 2A, or escape into the substrate. 
The energy carried by the escaping phonons is completely lost 
from the quasiparticle system. At low temperatures T ~ 0.177 
(where is the superconductor critical temperature), the 
downconversion process happens very fast (< 10 ns, as the 
quasiparticle-phonon scattering lifetime is much shorter 
than the recombination lifetime |^|^]) so it is the long-lived 
low energy quasiparticles of energy 7? rs A which primarily 


determine the detector response. We define an associated 
quasiparticle generation efficiency rj, given by the fraction 
of the absorbed photon energy that remains detectable as 
excess low energy quasiparticles. For very high energy photons 
(hi/ ^ Ud), rj = 0.6 is commonly used but over the 

moderate energy range (2A <hv < lOA - THz spectrum) of 
signal photons studied here, p varies significantly [^]. 

In [U we introduced a method for calculating the steady 
state, nonequilibrium quasiparticle and phonon distributions 
in superconducting thin films with simultaneous above-gap 
(signal) and sub-gap (readout or probe) photon illumination, 
by solving the nonlinear Chang & Scalapino kinetic equa¬ 
tions [p^. Using this method, the effect of uniform, constant 
absorption of sub-gap photons 0 0 and moderate energy 
above-gap photons [^, has been quantified. A key result of 
that work was an analytical relationship between the effective 
quasiparticle temperature T*n and absorbed power P (at a 
single photon frequency v). 


r]{v,P,Tb) (1 -I- n/Tpb) ^ 


( -2A{Tb) \ 
\ ksTb j 


( 1 ) 


Here the effective quasiparticle temperature Tf is defined as 
the temperature of the thermal distribution which has the same 
total number of quasiparticles as the steady-state nonequilib¬ 
rium distribution of interest. Using Tf in equations that as¬ 
sume thermal quasiparticle distributions is often sufficient for 
calculating key characteristics such as surface impedance [|]]. 
Equation (|^ therefore can be used to calculate the effective 
temperature from absorbed power (or vice versa) using only 
a material dependent constant Eg, derived from fitting the 
effective temperatures of the calculated nonequilibrium distri¬ 
butions; and the material independent but power, temperature 
and frequency dependent rj, also originally calculated from the 
nonequilibrium distributions, r/ is the phonon escape time into 
the substrate; Tpb is the phonon pair breaking time, which at 
low temperatures T Tc is equal to the characteristic phonon 
lifetime Tq [^; and Tb is the substrate or heat bath temperature. 
Our most recent work JI^ ] has calculated Eg for a range of 
common materials (Al, Mo, Ta, Nb, NbN), and also calculated 







26TH INTERNATIONAL SYMPOSIUM ON SPACE TERAHERTZ TECHNOLOGY, CAMBRIDGE, MA, 16-18 MARCH, 2015 


rj in the sub-gap and above-gap frequency regimes, at a range 
of temperatures. 

In this work, we compare (||) with the complete solutions 
to the Chang & Scalapino equations at the typical absorbed 
powers, signal and readout frequencies, temperatures, and 
for a range of commonly used low-Tc superconductors. We 
also compare the frequency dependence of our calculated 
quasiparticle generation efficiency p to recent measurements 
of Ta KID response [Q at THz frequencies. 

II. Results 



Fig. 1. Quasiparticle effective temperature against absorbed sub-gap 
power Pprobe for Al, Ta, Nb, Mo, and NbN, from full nonequilibrium 
calculation (markers) and analytical expression (lines). Calculated with hvp = 
16 ireV, T(, = 0.1 Tc and tiJtq = 1. 

In Fig. [|, we compare the quasiparticle effective temperature 
calculated from (|^ (lines) - for sub-gap readout frequency 
Vp and varying readout power Pprobe - to calculated from 
the full nonequilibrium distributions for Al, Ta, Nb, Mo, and 
NbN (markers), using the values for Eg tabulated in [|1- 
There is excellent agreement for all readout powers considered. 
The material-dependent constant Eg scales with the zero- 
temperature superconducting gap energy of the material, so 
for the same absorbed power, a greater effective temperature 
change is seen in Mo than NbN. 

In Fig. H the effective temperature calculated from (j^ 
(lines) is compared to n calculated from the nonequilibrium 
distributions (markers) when varying the substrate temperature 
Tfc. As the substrate temperature increases, the same absorbed 
power causes a smaller increase in the quasiparticle effective 
temperature T^, as temperature is a nonlinear function of total 
quasiparticle number and energy. Equation (|^ is in reasonable 
agreement with the full calculation until T;, ss 0.8 Tc, where 
ksT ^ A(r). 

Fig. I compares the effective temperatures as a function 
of absorbed sub-gap power for selected phonon escape time 
ratios. The effect of phonon trapping is correctly taken into 
account by (j^ as shown by the close agreement. Increasing 
the phonon escape time t; causes a greater increase in quasi¬ 
particle effective temperature for the same absorbed power. 



Fig. 2. Quasiparticle effective temperature difference — T;, against 
substrate temperature for selected absorbed sub-gap powers Pprobe, from 
full nonequilibrium calculation (markers) and analytical expression (lines). 
Calculated for Al, with hvp = 16 /reV and Tijr^ = 1. 



Fig. 3. Quasiparticle effective temperature against absorbed sub-gap 
power Pprobe for selected phonon escape time ratios r; /t^, from full nonequi¬ 
librium calculation (markers) and analytical expression (lines). Calculated for 
Al, with hvp = 16 /reV, and Tj = 0.1 Tc. 


as the probability for phonons to escape instead of breaking 
Cooper pairs is a function of Ti/Tpb. 

Fig. Q compares the effective temperatures as a function 
of absorbed power when the photons are sub-gap (dashed 
line, o markers - frequency hvp, absorbed power Pprobe) and 
above-gap frequency (solid line, -f markers - frequency Vg, 
absorbed power Pgignai), showing the analytical expression (j^ 
(lines) reproduces the quasiparticle effective temperature from 
the full calculation (markers). At the same absorbed power, all 
parameters of ([|) are identical between the sub-gap and above¬ 
gap cases except for 77. For direct pair breaking, the required 77 
is calculated from the full nonequilibrium distributions using 
a set of modihed Rothwarf-Taylor rate equations [^. In 
the sub-gap case 77 = 772 a, the fraction of phonons escaping 
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Fig. 4. Quasiparticle effective temperature against absorbed power 
P in the sub-gap (dashed line, o markers) and above-gap (solid line, + 
markers) cases, from full nonequilibrium calculation (markers) and analytical 
expression (lines). Calculated for Al, with hup = 16 /ieV, hus = lOA, and 

n = 0.1 Tc. 

the thin film which have energy 17 > 2A. rj is constant 
with absorbed above-gap power [^], while it decreases with 
absorbed sub-gap power Sin . For the same absorbed power, 
above-gap power is more efficiently converted into excess 
quasiparticles than sub-gap power - so the above-gap absorbed 
power results in a higher effective temperature than the sub¬ 
gap absorbed power. 



Fig. 5. Quasiparticle generation efficiency 77 as a function of signal frequency 
Us as calculated from full nonequilibrium distributions at two different phonon 
escape time ratios (dotted and .solid lines), compared to measured FTS spectral 
response of a Ta KID from BlJ] (dashed line). 

Neto et al. have recently measured the frequency depen¬ 
dence of the response of a Ta KID, on a SiN membrane to 
increase the thermal isolation (and hence t;), using a Fourier 
Transform spectrometer (FTS) ij^]. Their data analysis of 
the normalized response took account of antenna effects and 
band-defining filters, but does not remove the superconducting 
absorption efficiency, which cuts on after Hi's = 2A. The 


remaining oscillations of period ^ 0.5A are interpreted as 
due to standing waves and not intrinsic to the KID [|l4|]- Fig- 1 
shows their response measurements (dashed line) along with 
our calculations of the quasiparticle generation efficiency rj at 
two different trapping factors (dotted line and solid line). Our 
high trapping factor calculation of p (t;/t^ = 15, solid line) 
shows excellent agreement with the measurements. The high 
trapping factor is as expected for a device on a membrane 
and so thermally isolated from the substrate. A more detailed 
comparison is in progress. In our view, these experimental 
measurements confirm that rj is indeed frequency dependent 
as calculated by our model. Using the calculations of rj, we 
may be able to differentiate between different phonon trapping 
factors based on the device response, or choose the phonon 
trapping factor to achieve the desired response. 

III. Conclusion 

We have shown that we are able to calculate the quasipar¬ 
ticle effective temperature resulting from uniform absorbed 
power using a simple analytical expression (||) in agreement 
with the calculated steady state nonequilibrium distributions 
from a full solution of the nonlinear kinetic equations, for 
device operating parameters (materials, temperatures, powers, 
and photon frequencies) typical of KIDs and similar super¬ 
conducting pair breaking detectors. The effective temperature 
provides a good approximation to the average quasiparticle 
lifetimes and surface impedance of the superconductor [|]]. 
Therefore ([^) allows significantly simpler inclusion of nonequi¬ 
librium effects in higher level calculations and device models. 
For example, quasiparticle heating due to readout power in 
KIDs, which may lead to hysteresis [[^], can be included 
without calculating the detailed distributions. We also show 
the calculated frequency dependent response to above-gap 
frequency photons (particularly in the THz range), represented 
by the quasiparticle generation efficiency rj, is in agreement 
with recent experimental measurements. 

As a next step, we are exploring using (^ and our calcula¬ 
tions of Es to implement a complete electrothermal model of a 
KID, as described in [|^. This will allow exploration of device 
response, electrothermal feedback and hysteresis phenomena 
in detail. We also note that though presented in this work 
in the context of KIDs, (j^ is applicable to all illuminated 
superconducting thin films, and so is relevant to other devices, 
for example superconducting qubits []l7[-p^], resonator multi¬ 
plexers for Transition Edge Sensors (TESs) [ |20| |, and thin film 
parametric amplifiers pi]]. 
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